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The process of deeply virtual Compton scattering will be shortly introduced, and the latest results from 
measurements at the HERA ep-collider at DESY will be given. In particular, the cross section has 
been measured with increased statistics at the collider experiments HI and ZEUS, while HERMES for 
the first time reports measurements of the t-dependence of the beam— charge asymmetry on hydrogen 
and of a beam-charge asymmetry on deuterium. 



Hard exclusive processes such as the 
Deeply Virtual Compton Scattering (DVCS) 
process j*p — > 7p can be expressed in 
terms of Generalized Parton Distributions 
(GPDs) l ' 2 . GPDs give a generalized de- 
scription of the partonic structure of the nu- 
cleon, where the ordinary Parton Distribu- 
tion functions (PDFs) and the nucleon form 
factors turn out to be the kinematic limits 
and the moments of GPDs, respectively 3 . 
Besides the fact that the DVCS process ap- 
pears to provide the theoretically most di- 
rect access to GPDs, it is also investigated in 
the framework of diffractive processes. While 
the DVCS process is similar to the diffrac- 
tive electroproduction of vector mesons, it 
avoids the theoretical complications of need- 
ing further non-perturbative information due 
to the vector meson distribution amplitude 
in the final state. Furthermore, it is also 
unique among the hard scattering processes 
in that DVCS amplitudes (magnitude and 
phase) can be determined. This is possible 
through a measurement of the interference 
between the DVCS and Bethe-Heitler (BH) 
processes, in which the photon is radiated 
from a parton in the former and from the 
lepton in the latter process. These processes 
have an identical final state, i.e., they are in- 
distinguishable and thus the photon produc- 
tion amplitude r is given as the coherent sum 
of the amplitudes of the DVCS {tdvcs) and 



BH (tbh) processes. The cross section for 
the exclusive leptoproduction of photons is 
then given as (see Ref. 4 for full expression) 

d XB dQ>d\t\ d^ " ]TBH]2 + |TWCS|2 + L 

(1) 

where xb represents the Bjorken scaling 
variable, — Q 2 the virtual-photon four- 
momentum squared and t the square of the 
four-momentum transfer to the target. The 
azimuthal angle <p is defined as the angle be- 
tween the lepton scattering plane and the 
photon production plane. At leading twist, 
the BH-DVCS interference term / can be 
written as 

I = ±[c[ cos^ReM - Pi s[ sin^ImM], (2) 

where +(-) denotes a negatively (positively) 
charged lepton beam with longitudinal polar- 
ization Pi. The DVCS amplitude M is given 
by a linear combination of the nucleon form 
factors with the so-called Compton form fac- 
tors 5 , which are convolutions of the twist-2 
GPDs with the hard scattering amplitude. 

In the following, the HI and ZEUS 
cross section measurements, which access 
the squared DVCS amplitude \rr>vcs\ , and 
the HERMES azimuthal asymmetry mea- 
surements, which access the DVCS amplitude 
directly via the interference term /, will be 
presented. 
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1 Measurements of Cross Sections 
at HI and ZEUS 

The collider experiments HI and ZEUS have 
a similar experimental setup and thus employ 
a similar method for their analysis. Since 
the scattered proton escapes detection, the 
event selection is based on the detection of 
two electromagnetic clusters and at most one 
reconstructed track. With the incoming pro- 
ton defining the forward direction, the can- 
didate events are subdivided into two sam- 
ples where the positron candidate is either in 
the forward/central "large-Q 2 " region, or in 
the rear "small-Q 2 " region. A simulation of 
the BH contribution that describes the BH- 
dominated former sample can then be used to 
subtract the BH contribution from the latter 
sample, which contains BH events as well as 
DVCS events plus additional background. At 
leading twist, this directly leads to the DVCS 
cross section since the (^-dependent terms in 
the interference term (see Eq. 2) vanish as 
the measurement is integrating over the az- 
imuthal angle. 

The latest results on the 7*p cross sec- 
tion for the DVCS process, as obtained by 
HI, are shown in Fig. 1 as a function of 
the photon-nucleon invariant mass W for 
— t < 1 GeV 2 . These exceed the published 
results 8 by almost a factor of 4 in statis- 
tics. The NLO QCD calculation in Fig. 1 is 
based on a GPD parametrization 6 ' 7 . Since 
GPDs reduce to ordinary PDFs in a cer- 
tain kinematic limit, as mentioned above, the 
PDFs according to MRST2001 and CTEQ6 
are in turn chosen to serve as an input for 
the model calculations. The data agree with 
both models within the theoretical uncer- 
tainties, which are due to the up-to-now 
unmeasured i-dependence of the cross sec- 
tion. Since theoretical predictions are ab- 
sent as well, an exponential ansatz e~ b '* is 
assumed with 5 GeV -2 < fe < 9 GeV -2 
and b = b (l - 0.15 log(Q 2 /2)) GeV~ 2 . The 
range for the i-slope covers the measured 
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Figure 1. Preliminary HI results on the 7*p — ► jp 
cross section as a function of W for (Q 2 ) = 8 GeV 2 . 
The full error bars include the systematic error added 
in quadrature to the statistical error. The data are 
compared to a theoretical prediction 6,7 based on 
GPDs. The uncertainty in the predictions (shaded 
bands) is due to the unknown t— slope which is as- 
sumed to be between 5 GeV -2 (upper bound) and 
9 GeV -2 (lower bound). 



range for light vector meson production, i.e., 
it was chosen under the assumption that the 
production of real photons and light vec- 
tor mesons have a similar i-slope. Clearly, 
the foreseen direct measurement of the t- 
dependence will be extremely beneficial. 

Fig. 2 shows the preliminary HI data as 
a function of Q 2 and W together with one 
of the results from the GPD based model 
calculation, but this time at a fixed value 
of bo = 7 GeV -2 . In addition, the already 
published results from HI 8 and ZEUS 9 and 
the results from a color dipole model calcula- 
tion 10 are shown. Based on pp and j*p cross 
sections, the color dipole model predicts the 
1*P ~ * IP amplitude only for forward scat- 
tering, that is at t = t m i n . Thus, in order to 
calculate the integrated cross section, an ex- 
ponential ^-dependence with a fixed slope of 
7 GeV -2 is assumed as well. The differences 
between the results of the GPD based model 
and the color dipole model become apparent 
in the comparison of the cross section as a 
function of W, where the data show a pos- 
sible discrepancy between the HI and ZEUS 
results. However, a direct measurement of 
the ^-dependence is needed before the data 
can become conclusive. 
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Figure 2. Preliminary HI result together with the 
published results from HI 8 and ZEUS 9 on the 
1*P 7P cross section as a function of Q 2 and 
W. The data are compared to theoretical predictions 
based on a GPD 6 ' 7 and on a color dipole model 10 , 
using a t-slope of 7 GeV -2 . 



2 Measurements of Azimuthal 
Asymmetries at HERMES 

In contrast to the squared BH and DVCS am- 
plitudes, the interference term I (see Eq. 2) 
does depend on the sign of the beam charge. 
Therefore a measurement of a cross section 
asymmetry with respect to the beam charge 
can isolate the real part of the interference 
term, while the imaginary part can be iso- 
lated with a polarized lepton beam (P; 
0). Measurements of azimuthal asymmetries 
with respect to the beam spin, accessing the 
imaginary part of M via a sin <p modulation, 
have been reported on hydrogen ' as well 
as on deuterium and neon 13 . The extrac- 
tion of an azimuthal asymmetry with respect 
to the beam charge, accessing the real part 
of M via a cos cj> modulation, is described in 
the following using hydrogen and deuterium 
targets. 

At the fixed-target experiment HER- 
MES, events were selected if they contained 
exactly one photon and one charged track 
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Figure 3. Beam-charge asymmetry Ac for the hard 
electroproduction of photons off the deuteron as a 
function of the azimuthal angle <f> for the exclusive 
sample. 



identified as the scattered beam lepton. The 
important kinematic requirements imposed 
on the scattered lepton were Q 2 > 1 GeV 2 
and W > 3 GeV. The photon was identi- 
fied by detecting an energy deposition in the 
preshower scintillator and in the calorimeter 
without an associated charged track. The po- 
lar angle f? 7 * 7 between the virtual and the 
real photon is required to be between 5 and 
45 mrad. Since the recoiling proton was not 
detected, events were selected if the missing 
mass M x of the reaction ep — > e^X corre- 
sponds to the proton mass. Due to the fi- 
nite energy resolution the exclusive region is 
defined as —1.5 < M x < 1.7 GeV based on 
signal-to-background studies using a Monte- 
Carlo simulation. 

The beam-charge asymmetry as a func- 
tion of 4> is calculated as 



A c {4>) = 



N+{cp)-N-{$) 

N+i<t>) + N-(<t>y 



(3) 



where N + and N~ represent the single pho- 
ton yields normalized to the number of de- 
tected DIS events using the positron and elec- 
tron beam, respectively. It is shown in Fig. 3 
for the exclusive sample collected on an unpo- 
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Figure 4. The cos <fi amplitude of the beam-charge 
asymmetry on hydrogen and on unpolarized or spin- 
averaged polarized deuterium as a function of — t for 
the exclusive sample. The GPD model calculations 
use a factorized or a Regge— inspired t-dependence 
with or without a D-term contribution. 



larized or spin-averaged polarized deuterium 
target. The indicated fit to the asymme- 
try describes the data quite well and indeed 
yields the expected cos</> behavior. Fig. 4 
shows the cos <j) amplitudes on deuterium and 
hydrogen as a function of — i, which are de- 
rived from the fit to the beam-charge asym- 
metry in each —t bin. For both targets, 
the signal only becomes sizeable for larger 
values of —t. While their slightly different 
behavior at large — t values can be due to 
incoherent scattering on the neutron in the 
deuteron, effects from coherent scattering on 
the deuteron can be expected in the first — t 
bin but are not apparent there. Note that 
an intermediate result of the hydrogen anal- 
ysis 14 with a preliminary i-averaged value of 
0.11 ± 0.04 (stat.) ±0.03 (sys.) was derived 
at a mean — t value of 0.27 GeV 2 , it thus can 
be approximately compared to the result in 
the third —t bin. 

The theoretical calculations shown in 
Fig. 4 are carried out at the average kine- 
matics of every —t bin. They are based on 
a GPD model developed in Refs. 15 ' 16 . The 
data appear to favor the model with the non- 



factorized t-dependence and a vanishing con- 
tribution from the so-called D-term. It is 
apparent that measurements of the beam- 
charge asymmetry have a large predictive 
power given the fact that the electron sample 
used was quite small. 
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